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Abstract ' 

This paper describes the avian wing geometry 
(Seagull, Merganser, Teal and Owl) extracted from non- 
contact surface measurements using a three-dimensional 
laser scanner. The geometric quantities, including the 
camber line and thickness distribution of airfoil, wing 
planform, chord distribution, and twist distribution, are 
given in convenient analytical expressions. Thus, the avian 
wing surfaces can be generated and the wing kinematics can 
be simulated. The aerodynamic characteristics of avian 
airfoils in steady inviscid flows are briefly discussed. The 
avian wing kinematics is recovered from videos of three 
level-flying birds (Crane, Seagull and Goose) based on a 
two-jointed arm model. A flapping seagull wing in the 3D 
physical space is re-constructed from the extracted wing 
geometry and kinematics. 

1. Introduction 

Inspired by bird flight, early aviation researchers have 
studied avian wings as the basics of developing man-made 
flight vehicles. This methodology is clearly seen in the work 
of Lilienthal [1] and Magnan [2], This may be partially the 
reason why early aircraft designers like the Wright brothers 
tended to use thin airfoils by simply simulating bird wings. 
However, this situation was dramatically changed since 
thick airfoils (such as Gottingen and NACA airfoils) 
designed based on theoretical and experimental methods of 
aerodynamics achieved much higher lift-to-drag ratio at 
Reynolds numbers in airplane flight. Thus, study of avian 
wings becomes a marginalized topic that only interests a few 
avian biologists and zoologists. Nachtigall and Wieser [3] 
measured the airfoil sections of a pigeon’s wing. Oehme 
and Kitzler [4] measured the planform of 14 avian wings and 
gave an empirical formula for avian wing planforms. 

Recently, there is renewed interest in low-Reynolds- 
number flight and flapping flight in the aerospace 
community due to the need of developing micro-air-vehicles 
(MAVs). Hence, it is worthwhile to revisit the problem of 
the geometry and aerodynamics of avian wings. In this 
paper, we measure the surface geometry of several avian 
wings using a 3D laser scanning system. Based on these 
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measurements, we extract the basic geometrical properties of 
a wing such as the camber line, thickness distribution, 
planform and twist distribution, and generate 3D wing 
surfaces. The aerodynamic performance of the avian wing 
airfoils in steady inviscid flow is calculated in comparison 
with typical low-Reynolds-number airfoils. We explain how 
to recover the avian wing kinematics from videos of a level- 
flying bird. The present paper provides useful data for 
further biomimetic study of low-Reynolds-number wings 
and flapping wings for MAVs. 

2. 3D Laser Scanner on FARO Arm 

Figure 1 shows a FARO Arm (FARO Technologies. 
Inc.) to which a NVision’s 3D non-contact laser scanner is 
attached for wing surface measurements. The FARO Arm is 
a high accuracy hand-held mechanical device with an 
exchangeable probe, that is used to measure objects and 
features to create data of a surface. When the NVision 3D 
Scanner is attached and aligned to the arm, the capability of 
acquiring high-density point cloud data of a surface becomes 
available. On the FARO arm, the position of the scanner 
relative to a given coordinate system is known accurately. 
The accuracy of the surface data is within 0.041mm and data 
can be given in the coordinate system chosen. It is the 
fastest, smallest, and lightest hand-held non-contact 
scanning system available. 

Operating with ModelMaker software, the system 
works on the principle of laser stripe triangulation. A laser 
diode and stripe generator is used to project a laser line onto 
the object. The line is viewed at an angle by cameras so that 
height variations in the object can be seen as changes in the 
shape of the line. The resulting captured image of the stripe 
is a profile that contains the shape of the object. The 
software processes video data to capture surface shape in 
real time at over 23,000 points per second. The NVision 
Scanner uses digital camera synchronization to ensure 
precise measurements. It can scan a large variety of 
materials and colors including black, and work in almost any 
lighting conditions. ModelMaker is Windows-compatible 
software that outputs data in a variety of CAD formats. This 
system can generate millions of data points. For illustration. 
Figure 2 shows data cloud of the surface of a seagull wing 
obtained using this system. In this study, we only use a 
subset of data, that is, wing-cross-section data at selected 
spanwise locations. 
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3. Data Processing 

The upper and lower surface of an airfoil are 
expressed as addition and subtraction of the camber line and 
thickness distribution, z upper =z U) + z U) and 

z, ower = Z( C )-Z (I) , respectively. To extract the mean 

camber line from measurements, we use the Birnbaum- 
Glauert camber line [5] 

5!ll = hl2^ ri ( 1 -r 1 )ts n (2T 1 -ir, ( 1 ) 

c c n= , 

where TJ = x/c is the normalized chordwise coordinate and 
z (c)max ' s the maximum camber coordinate, and C is the 
local wing chord. The thickness distribution is given by [5] 
z ilL= z 11I ^j^ (2) 

c C n=l 

where z (t)max is the maximum thickness coordinate (the 


maximum thickness is 2z (cjmax ). For a given set of 

measured data of wing contour, a rotation and translation 
transformation is first applied in order that the geometrical 
angle-of-attack becomes zero and the leading edge of the 
wing section is located at the origin of the local coordinate 
system. Therefore, the local wing chord, twist angle, 
z tomax ’ z (t)max ’ relative position of the leading and trailing 
edges can be determined. Next, using least-squares 
estimation, we can obtain the coefficients S n and A n in 
Eqs. (1) and (2). These quantities are functions of the 
normalized span wise coordinate £ = 2y / b , where b/2 is the 
semi-span of a wing in a sense of the orthographic 
projection. Table 1 shows the averaged coefficients for the 
camber line and thickness distribution for the Seagull, 
Merganser, Teal and Owl wings. Details of how to extract 
these coefficients from measurements are discussed in the 
following sections. 

The chord can be expressed as 


c 

171 


_£o_ 

b/2 


. ^Ok( £ ) + Fcorr( £ ) ] ’ 


(3) 


where c 0 is the root chord of a wing. The function F 0K ( £ ) 
is a correlation given by Oehme and Kitzler [4] for avian 
wings, which is defined as F 0K (tq) = l for £e [0,0.5] 
and F ok ( q ) = 4^( 1 - £ ) for ge [0.5, 1 ] . The correction 
function for the deviation of an individual wing from 

F ok (£) is F co J$) = Y 4 E n (2; n+2 -2; s ), where the 

n=l 


coefficients E n are to be determined. Table 2 shows the 
coefficients for the planform of the Seagull, Merganser, Teal 
and Owl wings. The maximum camber line and thickness 
coordinates z (c)max and z U)max can be described by 


appropriate empirical functions of !q = 2y/b . Similarly, the 
relative position and kinematics of the 1/4-chord line of a 
wing to the fixed body coordinate system can be described 


by a dynamical system ( x c/4 , y c/4 , z c/4 )/( b/2 ) = f c/4 (t ), 
where t is time. As an approximate model, an avian wing 
can be described as a multiple jointed rigid arm system and 
its kinematics can be determined [6, 7]. In this paper, for 
simplicity, we adopt a two-jointed rigid arm system to 
describe the 1/4-chord line of an avian wing rather than 
exactly simulating the more complicated skeleton structure 
of an avian wing. The local twist angle of the airfoil section 
around the 1/4-chord line can be given by 
0 = fe( 2x c/4 / b,2y t , 4 / b, 2z c/4 /b) . When the geometric 
and kinematic parameters in the above relations are given, a 
flapping wing can be computationally generated and the 
wing kinematics can be simulated. 

4. Avian Wing Geometry 
4.1. Seagull 

Figure 3 shows a photograph of the Seagull wing used 
for this study. The coefficients S n and A n in Eqs. (1) and 
(2) for the camber line and thickness distribution are 
extracted from measurements of the Seagull wing. It is 
found that they do not show the systematic behavior as a 
function of the spanwise location, as shown in Fig. 4. 
Particularly, a considerable variation in A n exists. The 
averaged values of S n along the span in Eq. (1) are 
Sj =3.8735, S 2 =-0.807 and S 3 =0.771 . The averaged 
coefficients A n in Eq. (2) for the thickness distribution are 
Aj = -15.246 , A 2 = 26.482 , A 3 = -18.975 and 
A 4 =4.6232. Figure 5 shows the normalized camber line 
and thickness distribution for the Seagull wing generated by 
using the above averaged coefficients. These distributions 
exhibit the averaged airfoil of the Seagull wing over 
% = 2y/b = 0.166 -0.772. For 2y/b> 0.772, the 
primaries are separated such that no single, continuous 
airfoil exists. The least-squares estimation residuals in 
fitting local airfoils z (c) /c and z (t) /c at different 
spanwise locations are shown in Fig. 6(a). Similarly, the 
deviations of the averaged z (c) /z (c)max and z ft) /z (t)max 

from the local profiles at different spanwise locations are 
shown in Fig. 6(b). 

As shown in Fig. 7, the maximum camber and 
thickness coordinates z lc)max and z (l)max are functions of the 

spanwise location % = 2y/b, which are empirically 
expressed as z (c)max / c = 0.14/( 1 + 1.333 g !A ) and 
Z( t )mu / c =0.1/( 1 + 3.546 £ ,A ) . Figure 8 shows the 

planform of the Seagull wing. The distribution of the wing 
chord, as shown in Fig. 9, can be described by Eq. (3) that is 
the Oehme and Kitzler’ s correlation F 0K (lq) plus a 

5 

correction function F corr { !q ) = ^ E n <^ n+2 -£* ) for local 

n=l 
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variation, where E, =26.08, E 2 =— 209.92, £ ? =637.21 , 
E 4 = -945.68 and E s = 695.03 . The ratio between the root 
chord and semi-span is c 0 /( b/ 2 ) = 0.388 . Figure 10 
shows the wing twist as a function of the spanwise location, 
which is expressed as an expansion of the Chebyshev 
polynomials ( T, = % , T 2 = 4% 3 - 35, , and 

3 

T 3 = 16% 5 - 20% 3 + 5 % ), i.e., twist (deg) = D n T„ ( % ) , 

n—1 

where D, =5.2788, D 2 =-4.1069 and D 3 = -1.8684. 
Here the positive sign of the twist denotes that the wing 
rotates against the incoming flow. Note that the wing twist 
presented here is not necessarily intrinsic because not only 
the twist may be changed in preparing the wing specimen, 
but also the twist is really a time-dependent variable during 
flapping. Using the above relations obtained from 
measurements, we generate the surface of the Seagull wing 
shown in Fig. 1 1 , where a simple two-jointed arm model is 
used for the 1/4-chord line. Also, we assume that the airfoil 
section remains the same near the wing tip while the 
maximum thickness decreases even though the real wing has 
separated primaries near the wing tip. 

4.2. Merganser 

Figure 12 shows a photograph of the Merganser wing 
used for this study. The coefficients S n and A n for the 
camber line and thickness distribution of the Merganser 
wing are shown in Fig. 13. The averaged values of S n 
along the span are S,= 3.9385, S, =0.7466 and 
S 3 = 1.840. The averaged coefficients A n for the thickness 
distribution are Aj= -23.1743, A 2 = 58.3057 , 
A ; = -64.3674 and A 4 — 25.7629. Figure 14 shows the 
normalized camber line and thickness distribution for the 
Merganser wing generated by using the above averaged 
coefficients. The wing thickness is very small (considered 
to be zero) near the trailing edge (x/c > 0.9). The least- 
squares estimation residuals in fitting local airfoils z (c) /c 

and z (t) /c at different spanwise locations are shown in Fig. 
15(a). The deviations of the averaged z (c) 7 z (c)max and 
z ut / Z U) max fr° m die local profiles at different spanwise 

locations are shown in Fig. 15(b). 

Figure 16 shows the maximum camber and thickness 
coordinates z (c)max and z U)ma as a function of the spanwise 

location % = 2y/b along with the empirical expressions 
z (c)nax /c = 0.14 /( 1+1.333 %' 4 ) and 

z u)max / c = 0.05 /( 1 + 4 %' 4 ). Figure 17 shows the 

planform of the Merganser wing. The distribution of the 
wing chord is shown in Fig. 18 along with the results given 
by Eq. (3) where the coefficients in F con _ ( % ) are E, =39.1, 


E 2 = -323.8 , Ej = 978. 7, E 4 =-141 7.0 and 
E 5 =1001.0. The ratio between the root chord and semi- 
span is c 0 /(b/ 2 ) = 0.423 . Figure 19 shows the wing twist 
as a function of the spanwise location, which is expressed as 
an expansion of the Chebyshev polynomials ( T, = % , 

T 2 = 4% 3 - 3 % , and T 3 =16% 5 -20% 3 +5%), i.e., 

3 

twist (deg) = ^ D n T n ( % ) , where D, = 30. 9953 , 

n-1 

D 2 = -3.2438 and D 3 = -0.2076 . Figure 20 shows the 
surface of the Merganser wing generated using the above 
relations. 

4.3. Teal 

Figure 21 shows a photograph of the Teal wing used 
for this study. The averaged values of S n along the span for 
the camber line are S , =3.9917 , S 2 =-0.3677 and 
S 3 = 0.0239. The averaged coefficients A n for the 
thickness distribution are = i./ou4, A 2 = — 13.68/5 , 
A } = 18.276 and A 4 =-8.279. Figure 22 shows the 
normalized camber line and thickness distribution for the 
Teal wing generated by using the above averaged 
coefficients. The least-squares estimation residual in fitting 
local airfoils z (c) /c and z (l) /c is less than 0.003. The 

deviations of the averaged z (c) /z (c)max and z u) /z u)ma 

from the local profiles are less than 0.1 and 0.2, respectively. 

Figure 23 shows the maximum camber and thickness 
coordinates z (c)max and z lt)max as a function of the spanwise 

location % = 2y/b along with the empirical expressions 
z (c)max /c = 0.11/(l + 4% 14 ) and 

z (t>mnx = 0.057(1 + 4 % 14 ) . Figure 24 shows the 

planform of the Teal wing. The distribution of the wing 
chord is shown in Fig. 25 along with the results given by Eq. 
(3) where the coefficients in F con (%) are E,=-66.1, 
E 2 = 435.6, E 3 =-1203, E 4 =1664.1 and E s =-1130.2. 
The ratio between the root chord and semi-span is 
c 0 /(b/ 2 ) = 0.545 . The wing twist is less than 2 degrees 
along the span. Figure 26 shows the surface of the Teal 
wing generated using the above relations. 

4.4. Owl 

Figure 27 shows a photograph of the Owl wing. The 
averaged values of S n along the span for the camber line are 
S, =3.9733, S 2 =-0.8497 and S 3 =-2.723. The 
averaged coefficients A n for the thickness distribution are 
A, =-47.683, A 2 = 124.5329, A 3 =-127.0874 and 
A 4 = 45.876 . Figure 28 shows the normalized camber line 
and thickness distribution for the Owl wing generated by 
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using the above averaged coefficients. The least-squares 
estimation residual in fitting local airfoils z u) / c and 

z (l) /c is less than 0.006. The deviations of the averaged 
Z (c) /Z (c)max and z (l) /z (nmax from the local profiles are 
less than 0.1 and 0.2, respectively. Interestingly, the Owl 
wing is very thin over x / c = 0.3— 1.0 (it is a single layer of 
the primary feathers) and the thickness distribution is mainly 
concentrated in the front portion of the airfoil. The wing 
thickness is considered to be zero near the trailing edge ( x/c 
>0.9). 

Figure 29 shows the maximum camber and thickness 
coordinates z (c)max and z„ )max as a function of the spanwise 

location £ = 2y/b along with the empirical expressions 
z (c)max / c = 0.04[ l+tanh( 1.8g-0.5 )] and 

z ( , >ma x / c - 0.04 /( 1 + 1.78% 14 ). In contrast to other wings 

described before, the maximum chamber coordinate for the 
Owl wing increase along the span. Figure 30 shows the 
planform of the Owl wing. The distribution of the wing 
chord is shown in Fig. 31 along with the results given by Eq. 
(3) where the coefficients in F corr (t;) are E, =6.3421, 
E 2 = -7.51 78, Ej = -70. 9649 , E 4 = 188.0651 and 
E s - -160.1678 . The ratio between the root chord and 
semi-span is c 0 /( b / 2 ) -0.677 . The wing twist is less 
than 2 degrees along the span. Figure 32 shows the surface 
of the Owl wing generated using the above relations. 

5. Aerodynamic Characteristics of Avian Airfoils in 
Steady Inviscid Flows 

Figure 33 shows typical wing sections of the Seagull, 
Merganser, Teal and Owl at 2y / b-0.4 . These airfoils are 
highly cambered. The inviscid pressure coefficient C p 

distributions at four angles of attack are shown in Figs. 34, 
35, 36 and 37. Figure 38 shows the sectional lift coefficient 
based on unit chord as a function the angle of attack (AoA) 
for the Seagull, Merganser, Teal and Owl wings. These 
results are obtained by using the inviscid/viscous flow 
analysis code XFOIL for airfoil design [8], which roughly 
indicate the aerodynamic characteristics of these airfoils. 
The pressure distributions on the upper surfaces of the 
Seagull and Merganser wings are relatively flat when AoA is 
less than 5 degrees. The sectional lift coefficients of at zero 
AoA for both are larger than one. Figure 39 shows the 
sectional lift coefficient distributions along the wing span for 
these wings at AoA = 0 degree. Based on the sectional life 
coefficient c, , we can estimate the normalized circulation 
distribution r( y ) / T 0 - [ c( y )/ c g ][c,( y )/ c w ] shown 
in Fig. 40, where the subscript ‘0’ denotes the value at the 
wing root (2y / b = 0 ). The Seagull and Merganser wings 
have the almost same normalized circulation distributions. 
The Owl airfoil is particularly interesting, that is basically a 


thin wing with a thickness distribution concentrated mainly 
near the leading edge. Unlike other wings, the c, 
distribution for the Owl wing has an increasing behavior as 
the wing span because the maximum camber coordinate 
Z( C )max increases. As a result, the normalized circulation 

distribution has a special shape as indicated in Fig. 40. We 
do not know whether the thin Owl wing and the associated 
aerodynamic properties are related to quiet flight of an owl 
[9]. Clearly, the aerodynamic and aeroacoustic implications 
of the thin Owl wing are worthwhile to be investigated 
further. 

The Seagull and Merganser airfoils are similar to the 
high-lift low Reynolds number airfoil SI 223 described by 
Selig et al. [10]. Figure 41 shows the S1223 airfoil along 
with the Seagull and Merganser airfoils with the same the 
maximum camber line and thickness coordinates 
(z (c)max /c = 0.0852 and z (l)max /c = 0.0579). Figures 42 

and 43 shows a comparison of the pressure coefficient 
distributions between the SI 223, Seagull and Merganser 
airfoils. These pressure distributions are similar, but the 
SI 223 airfoil has lower pressure on the upper surface near 
x/c = 0.2 and trailing edge. The sectional lift coefficient as a 
function of AoA for these airfoils is shown in Fig. 44. When 
AoA increases beyond a certain value (about 10 degrees), 
laminar flow separation will take place near the leading edge 
in a Reynolds number range for birds (4x1 0 4 to 7xl0 5 ) 
[11,12]. The separated flow may be reattached due to 
transition to turbulence that can be facilitated by using 
artificial boundary layer tripping. Detailed calculation of the 
separated/reattached flow on these airfoils requires a Navier- 
Stokes (N-S) solver with accurate transition and turbulence 
models; computation based on the N-S equation especially 
for the unsteady flow field around a flapping wing is a topic 
in further study. Here, we do not intend to conduct such 
computation without reliable experimental data for 
comparison. Nevertheless, experimental data for the SI 223 
airfoil [10] provide a good reference (in a qualitative sense) 
for the behavior of the Seagull and Merganser airfoils at 
high angles of attack. 

6. Avian Wing Kinematics 

6.1. Front-Projected 1/4-Chord Line 

For simplicity, we consider the kinematics of a 
flapping wing as a superposition of the motion of the 1 74- 
chord line of the wing and relative rotation of local airfoil 
sections around the 1/4-chord line. From videos of a level- 
flying bird taken by a camera viewing directly the front of 
the bird, we are able to approximately recover the front- 
projected profiles of the 1 /4-chord line of the wing at a 
sequence of times. Figure 45 shows a typical image of a 
level-flying crane viewed directly from the front and a local 
coordinate system used for describing the profiles. A front- 
projected wing in images is a line with a finite thickness that 
is approximately considered as the front-projected 1 /4-chord 
line. The profile of the front-projected 1/4-chord line of a 
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flapping wing can be reasonably described by a second- 
order polynomial 


b/2 


= A,( cot ) 


! 3u± 
b/2 


+ A 2 (6)t ) 


r i ^ 

b/2 


(4) 


where the coefficients and the semi-span are given by the 
Fourier series as a function of the non-dimensional time a t 


( co is the circular frequency of flapping) 

2 

A t (a)t ) = C I0 +'Y J [c, n sinfneot )+B In cos( ncot )], 

n=l 

2 

A 2 ( cot ) = C 20 +£[C 2 „ sin( ncot )+B 2n cos( ncot )], 

n=l 

b( cot ) / 2 r 1 

/I = sm <n<»t)+B bn cos(ncot) \. 

fticix( b/2 ) n _j 

(5) 

Here, b/2 is defined as the semi-span of an orthographically 
projected flapping wing on the horizontal plane. Therefore, 
b/2 is a time-dependent function in a flapping cycle. The 
maximum value of b/2 is achieved roughly at the moment 
when a flapping wing is parallel to the horizontal plane. We 
assume that cot =0 corresponds to the position of a wing at 
the beginning of the down-stroke (or the end of the up- 
stroke) (see Fig. 45). 


6.1.1. Crane 

A time sequence of images of a level-flying crane 
taken by a camera directly from the front of the bird are 
obtained by digitizing a clip of the video “The Life of Birds” 
produced by BBC. The profiles of the front-projected wing 
(or 1/4-chord line) are obtained by manually tracing the 
wing in digitized images. Eq. (4) is used to fit data of the 
successive profiles and the coefficients in Eq. (5) are 
determined. Figure 46 shows the measured profiles of the 
front-projected 1/4-chord line of a flapping wing of a flying 
crane and the corresponding polynomial fits at six instants 
(an interval of 2k/ 5) in a flapping cycle cots [ 0/2 k ] . 
The profiles can be reasonably described by a second-order 
polynomial Eq. (4) with the time-dependent coefficients. 
Figures 47 and 48 show data of the coefficients in Eq. (5) 
and the orthographically projected semi-span b/2 that are fit 
by the Fourier series, respectively. The coefficients in Eq. 
(5) extracted from measurements for a flapping wing of a 
crane are 

C l0 = 0.3639 , C„ = -0.2938 , B,, = 0.4050 , 

C, 2 = -0.0465 , B I2 = -0.0331 ; 

C 20 = -0.4294 , C 2I = 0.4469 , B 2I = 0. 1442 , 
C 22 = 0.0135 , B 22 = 0.0691 ; 

C h0 =0.839, C bI =0.0885, B bI =0.0301, C h2 =-0.0888, 
B b2 =-0.0407 . 

Figure 48 shows that the orthographically projected 
semi-span b/2 on the horizontal plane varies with time. At 


cot =0 , the position of the wing is at the beginning of the 
down-stroke. The wing is approximately parallel to the 
horizontal plane at cot = 2 and b/2 reaches the maximal 
value. The minimal value of b/2 is at cot = 3.9 . The down- 
stroke spans about 62% of a flapping cycle while the up- 
stroke takes 38% of a cycle. The variation of b/2 with time 
depends on not only the orthographic projection, but also a 
change of the wing planform due to wing extension and 
folding during flapping. We calculate the arc length of the 
front-projected 1 /4-chord line as a function of time by using 
Eqs. (4) and (5). In fact, a change in the arc length of the 
front-projected 1 /4-chord line represents a change of the 
wing planform due to wing extension and folding. Figure 49 
shows the arc length of the projected 1 /4-chord line as a 
function of time for the flapping crane, seagull and goose 
wings. For a crane, its wing is most extended at cot =2.1 
while it is most folded at cot = 4. The normalized arc 
length of the front-projected 1 /4-chord line is described by 
the Fourier series 

¥ / \ ? 

is • | l If l f ^ , r 

— = C L0 + Z_\ C Ln sin( ncot 1+B^, cos( ncot )]. 

W L pmj ) 

(6) 

For the flapping crane wing, the coefficients in Eq. (6) are 
C L0 = 0.9310 , C L1 = 0.0359 , B u = 0.0111 , 

C L2 =-0.0675, B u =-0.0093. 

This result will be used later to re-construct the wing 
kinematics based on a two-jointed arm model. 

6.1.2. Seagull 

Similarly, a time sequence of images of a flying 
seagull (acquired from Oceanfootage.com) is processed and 
the profiles of the front-projected 1/4-chord line are 
recovered. The coefficients in Eq. (5) extracted from 
measurements for a flapping wing of a seagull are 
C w = 0.3756 , C„ = -0.3242 , B„ =0.1920, 

C I2 =0.0412, B 12 =-0.1095 ; 

C 20 = -0.4674 , C 2I = 0.3631 , B 2I = 0.2884 , 

C 22 = -0.0661 , B 22 =0.0553 ; 

C M =0.7978, C bI =0.1751, B bl =0.0461, C b2 =0.0042, 
B b2 =-0.0218. 

Figure 50 shows the measured profiles of the front-projected 
1/4-chord line of a flapping wing of a flying seagull and the 
corresponding polynomial fits at six instants (an interval of 
2k/ 9) in a flapping cycle. Figures 5 1 and 52 show data of 
the coefficients in Eq. (5) and the orthographically projected 
semi-span b/2 that are fit by the Fourier series, respectively. 

For the normalized arc length of the front-projected 
1/4-chord line of the flapping seagull wing, the coefficients 
in Eq. (6) are 

C L0 = 0.8718 , C LI = 0.1420 , B u = -0.0111 , 

C L2 =0.0190. B u =0.0113. 
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As shown in Fig. 49, the flapping seagull wing is most 
extended at cot = 1.3 while it is most folded at cot ~ 5 . 

6.1.3. Goose 

A time sequence of images of a flying bar-headed 
goose from the documentary “Winged Migration” is 
processed and the profiles of the front-projected 1/4-chord 
line are recovered. The coefficients in Eq. (5) extracted 
from measurements for a flapping wing of a level-flying 
goose are 

C, 0 = 0.4511 , C„ = -0.2819 , B„ = 0.3008 , 

C , 2 = 0.0254, B, 2 = -0.0835 ; 

C 20 = -0.4605 , C 2I = 0.4516 , B 21 = 0.1912 , 

C 22 =-0.0845, B 22 =0.1154; 

C b0 = 0.8999 , C bI = 0.0666 , B h , = 0.0126 , 

C b2 = -0.0505 , B b2 = -0.0095 . 

Figure 53 shows the measured profiles of the front-projected 
1 /4-chord line of a flapping wing of a flying goose and the 
corresponding polynomial fits at six instants (an interval of 
n / 5 ) in a flapping cycle. Figures 54 and 55 show data of 
the coefficients in Eq. (5) and the orthographically semi- 
span b/2 that are fit by the Fourier series, respectively. 

For the normalized arc length of the front-projected 
1/4-chord line of the flapping goose wing, the coefficients in 
Eq. (6) are 

C L0 = 0.9948 , C LI = 0.0013 , B u = -0.0013 , 

C L2 =-0.0083, B u =0.0122. 

As shown Fig. 49, the normalized arc length of the front- 
projected 1/4-chord line of the flapping goose wing does not 
vary much compared with the flapping crane and seagull 
wings. This means that relatively speaking the goose wing 
does not extend and fold much during flapping. 

6.2. Two- Jointed Arm Model 

In general, the skeleton structure is described as a 
three-jointed arm system. Figure 56 is an X-ray image 
showing the skeleton structure of a seagull wing. However, 
for level flapping flight, the wing kinematics can be 
simplified. In this case, to describe the 1/4-chord line of a 
flapping wing, we use a two-jointed arm model that consists 
of two rigid jointed rods. As shown in Fig. 57, Rod 1 rotates 
around the point O, in a body coordinate system where the 
origin O, is located at the wing root and the plane YO,Z is 
defined as the rotational plane of Rod 1. Thus, the motion of 
Rod 1 has only one degree of freedom and the position of 
Rod 1 is given by the flapping angle yf, . In contrast, the 
motion of Rod 2 has two degrees of freedom, which is given 
by the angles yr 2 and <j> 2 . In Fig. 57, the line 0 2 T' is the 
orthographic projection of the Rod 2 (or the line 0 2 T ) on 
the plane YO,Z . The angle ys 2 is the angle between Rod 1 
and the line 0 2 T' on the plane YO,Z , which basically 


determines the flapping magnitude of Rod 2 relative to Rod 
1. The angle <j) 2 is the angle between Rod 2 and the line 
0 2 T' , which describes the extension and folding of a wing 
(the outer portion of a wing). Figure 58 shows the projected 
views of a two-jointed arm system. In Fig. 58(c), the angle 
<p 2± = <p 2 / cos( y/, - iff 2 ) is the orthographic projection of 
the angle <p 2 on the horizontal plane XO, Y . The simple 
two-jointed arm model allows the recovery of 3D kinematics 
of a flapping wing from measurements of the front-projected 
1/4-chord line. In addition, it is a straightforward model for 
designing a mechanical flapping wing. 

The coordinates of the end point 0 2 of Rod 1 are 


x o 2 = 0 • y o 2 = L i cos ( Vi )> z o 2 = L , sin < ¥, ) > (7) 

where L , is the length of Rod 1 . The position of Rod 1 is 
described by 

X =0 

„ „ , ( 8 ) 

Z = Y tan( \j/ 1 ) 

where Y e [0,L, cos( yf,)]- The position of Rod 2 is given 
by 

✓ \ 

02 


-(Y-Y 02 )tan 


cos( y/,-y/ 2 ) 


(9) 


Z = Z 02 +(Y- Y 02 )tan {yf, - yf 2 ) 
where Y e [ L, cos( yf , ),b / 2 ] . Note that b/2 is the 

orthographically projected semi-span on the horizontal plane 
XOjY . Therefore, we know that the projected semi-span is 
b/2 = L, cos(yf, )+L 2 cos(tp 2 )cos(yf,-yf 2 ). In a two- 
jointed arm system, the normalized arc length of the front- 
projected 1/4-chord line is 


L p roj(<»0 . 
max( L proj ) 


where r, = L, / max( L proj 
the relative lengths of Rod 1 and Rod 2 


: r,+r 2 cos(<p 2 ), (10) 

) and r 2 =L 2 / max( L proj ) are 


6.3. Recovery of the Angles yf , , yf 2 and (p 2 

A two-jointed arm model uses two pieces of straight 
line to approximate the profile of the 1 /4-chord line of a 
wing. Since the flapping angles yf, and yf 2 are on the plane 
YO,Z , they can be estimated directly from the measured 
profile of the front-projected 1/4-chord line, Eq (4), when 
r,=L,/ max( L pwj ) and r 2 = L 2 / max( L proj ) are given. 

The angle <j) 2 can be extracted from the measured arc length 
of the front-projected 1/4-chord line using Eqs. (10) and (6). 
Figures 59, 60 and 61 shows the recovered angles yf , , yf 2 
and (f> 2 as a function of time for the flapping crane, seagull 
and goose wings, respectively. 

The angles yf, , yf 2 and (j) 2 are expressed as the 
Fourier series 
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y/,(( 0 t) = C v]0 +'Y^C v , n sin(na)t)+B ¥ln cos(nm)], 

n~l 

V 2 ((0t) = C v2O +'Y J [c ¥ 2 n sin( rnot)+B ¥2n cos(na)t)], 

n=I 

<p 2 (cot)= C, 20 + \c, 2n sin( nan )+B 02n cos( neot )]. 

n=l 

(H) 

The estimated coefficients in Eq. (11) for the crane, seagull 
and goose wings are given below. Here, we assume that 


r ; = 0.5 and r 2 = 0.5 . The units of the angles y/j , y/ 2 and 
<j> 2 in Eq. (1 1 ) are in degrees. 

Crane 

C ¥l0 =8.3065, C ¥ „ =-4.4519, 

C vI2 =-1.8092 , B ¥l2 =-0.5889 ; 

B ¥ll =25.3910, 

C ¥20 = 17.0661, C ¥2I =-17.3404, 

C ¥ 22 =-3.7029, B„ 22 = —4.5122 ; 

B ¥2I =-4.0066 , 

C 02O =32.231 1 , C, 2I =-8.6004, 

C 022 =15.2213, B t22 =2.6910. 

B 02I =-0.34280, 

Seasull 

C ¥l0 = 8.4654 , C ¥lI = -8.5368 , 

C ¥l2 = 1.0898 , B ¥l2 = -4.5880 ; 

B ¥ „ =17.8798, 

C ¥20 = 17.3083 , C ¥2I = -11.0122 , 
C ¥22 =1.3128, B ¥22 =-3.0183 ; 

B ¥21 =-9.6131, 

C 02O = 38.41 79 , C, 2I = -28.0553 , 

C $22 =—41032 , B t22 =3.0125. 

B 02 , =0.7664, 

Goose 

C ¥I0 = 12.2528 , C ¥lI = -3. 7150 , 

C ¥l2 =-0.6432 , B ¥l2 =-2.3054 ; 

B ¥U = 21.1873 , 

C ¥20 = 20.0863 , C ¥ 2j = -18.6807 , 

C ¥22 =1.3467, B ¥22 =-6.1507 ; 

B ¥2 , =-7.3848, 

C, 20 = 13.5235 , C, 2I = -0. 7494 , 

C, 22 =4.3138 ', B,22 =-6.3023. 

B, 2 , =1.2524, 


6.4. Reconstruction of a Flapping Wing 

After the wing geometry (the airfoil section, 
planform, and twist distribution) and the kinematics of the 
1/4-chord line of a wing are given, a flapping wing can be 
re-constructed in the 3D physical space by superimposing 
the airfoil sections on the moving 1 /4-chord line. Note that 
the wing twist distribution in flapping is not recovered in 
this paper. Measurements of the dynamical wing twist 
distribution require considerable videogrammetric 
processing on a time sequence of images taken from two 


cameras simultaneously viewing a flapping wing of a level- 
flight bird on which a sufficient number of suitably 
distributed targets are attached. In computational 
simulations, the wing twist can be treated as a variable to 
achieve the maximum aerodynamic efficiency. Here, we 
simply assume that the wing twist is fixed during flapping. 
Using Eqs. (1), (2), (3), (8), (9) and (11) with the known 
coefficients for a seagull wing, we re-construct a flapping 
seagull wing at different instants as shown in Fig. 62. 

6. Conclusions 

Using a 3D laser scanner, we have measured the 
surface geometry of the Seagull, Merganser, Teal and Owl 
wings. From measurements, the airfoil camber line, airfoil 
thickness distribution, wing planform and twist distribution 
are extracted. The accuracy of metric measurements using 
the laser scanner is about 0.041mm. The residual of least- 
squares fitting for an airfoil section is about 2-10xl0 3 in 
terms of the normalized coordinate z/c. The estimated 
coefficients for the camber line and thickness distribution do 
not exhibit a systematic behavior along the wing span. 
Thus, the averaged values of these coefficients along the 
wing span are given, which define the averaged airfoil for an 
avian wing. The deviation of the local airfoil camber line 
and thickness distribution from the averaged ones is about 5- 
20% of their maximum value. The Seagull and Merganser 
airfoils are similar to high-lift low Reynolds number airfoils. 
The Teal airfoil has a relatively symmetric thickness 
distribution around the mid-chord. The Owl airfoil is very 
thin over 0.3 -1.0 chord and the thickness distribution is 
mainly concentrated in the front portion of the airfoil. 
Unlike other wings, the Owl wing has a special circulation 
distribution along the wing span. 

We consider the kinematics of a flapping wing as a 
superposition of the motion of the 1/4-chord line of the wing 
and relative rotation of local airfoil sections around the 1/4- 
chord line. The profiles of the front-projected 1/4-chord line 
at different instants are measured from videos of a level- 
flying bird. Then, based on a two-jointed arm model, the 
kinematics of the 1/4-chord fine in the 3D physical space is 
recovered for the flapping Crane, Seagull and Goose wings. 
The relevant quantities of the wing kinematics are given in 
convenient analytical expressions. The wing geometry and 
kinematics given in this paper are useful for the design of 
flapping MAVs and experimental and computational studies 
to understand the fundamental aerodynamic aspects of 
flapping flight. 
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Table 1. The Coefficients for Avian Airfoil 



Seagull 

Merganser 

Teal 

Owl 

S, 

3.8735 

3.9385 

3.9917 

3.9733 

S 2 

-0.807 

0.7466 

-0.3677 

-0.8497 

S., 

0.771 

1.840 

0.0239 

-2.723 

A] 

-15.246 

-23.1743 

1.7804 

-47.683 

A 2 

26.482 

58.3057 

-13.6875 

124.5329 

A3 

-18.975 

-64.3674 

18.276 

-127.0874 

a 4 

4.6232 

25.7629 

-8.279 

45.876 


Table 2. The Coefficients for Wing Planform 



Seagull 

Merganser 

Teal 

Owl 

E, 

26.08 

39.1 

-66.1 

6.3421 

E, 

-209.92 

-323.8 

435.6 

-7.5178 

Ei 

637.21 

978.7 

-1203.0 

-70.9649 

e 4 

-945.68 

-1417.0 

1664.1 

188.0651 

e 5 

695.03 

1001.0 

-1130.2 

-160.1678 



Figure 1. 3D laser scanner and FARO arm for wing surface 
measurements. 




Figure 3. The Seagull wing. 
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Figure 4. (a) The coefficients for the camber line, (b) The 
coefficients for the thickness distribution for the Seagull 
wing. 
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Figure 5. The camber line and thickness distribution of the 
Seagull wing. 
x10 3 



Spanwise Position Normalized by b/2 

(a) 



(b) 

Figure 6. (a) Least-squares residuals of fitting the airfoil 
sections, (b) Deviation of local profiles from the averaged 
profile for the Seagull wing. 


o 



Figure 7. The maximum camber and thickness coordinates 
as a function of the spanwise location for the Seagull wing. 



Figure 8. The planform of the Seagull wing. 
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Figure 9. The chord distribution of the Seagull wing. 
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Figure 10. The twist distribution of the Seagull wing. 




Figure 12. The Merganser wing. 
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Figure 13. (a) The coefficients for the camber line, (b) The 
coefficients for the thickness distribution for the Merganser 
wing. 


Figure 1 1 . The generated surface of the Seagull wing. 
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Figure 14. The camber line and thickness distribution of the 
Merganser wing. 
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Figure 15. (a) Least-squares residuals of fitting the airfoil 
sections, (b) Deviation of local profiles from the averaged 
profile for the Merganser wing. 
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Figure 16. The maximum camber and thickness coordinates 
as a function of the spanwise location for the Merganser 
wing. 
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Figure 17. The planform of the Merganser wing. 
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Figure 18. The chord distribution of the Merganser wing 
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Figure 22. The camber line and thickness distribution of the 
Teal wins. 


Normalized Spanwise Coordinate. 2y/b 
Figure 19. The twist distribution of the Merganser wing 
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Figure 23. The maximum camber and thickness coordinates 
as a function of the spanwise location for the Teal wing. 


Figure 20. The generated surface of the Merganser wing 
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Normalized Spanwise Coordinate, 2y/b 


Figure 24. The planform of the Teal wing. 



Normalized Spanwise Coordinate, 2y/b 
Figure 25. The chord distribution of the Teal wing. 



Figure 26. The generated surface of the Teal wing. 
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Figure 28. The camber line and thickness distribution of the 
Owl wing. 
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Figure 29. The maximum camber and thickness coordinates 
as a function of the spanwise location for the Owl wing. 
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Figure 30. The planform of the Owl wing. 
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Figure 33. Airfoil sections of the avian wings at 2y/b = 0.4. 
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Figure 35. The pressure coefficient distributions of the 
Merganser wing at different angles of attack. 
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Figure 34. The pressure coefficient distributions of the 
Seagull wing at different angles of attack. 


Figure 36. The pressure coefficient distributions of the Teal 
wing at different angles of attack. 
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Figure 37. The pressure coefficient distributions of the Owl 
wing at different angles of attack. 


Figure 39. The sectional lift coefficient distributions along 
the wing span for the Seagull, Merganser, Teal and Owl 
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Figure 38. The sectional lift coefficient as a function of the 
angles of attack. 


Figure 40. The normalized circulation distributions along the 
wing span for the Seagull, Merganser, Teal and Owl wings. 
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Figure 45. An image of a level-flying crane viewed by a 
camera directly from the front. 
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Figure 46. The profiles of the front-projected 1/4-chord line 
of the flapping crane wing at different instants. 
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Figure 48. The orthographicaliv projected semi-span b/2 
normalized by max(fc/2) for the flapping crane wing as a 
function of time. 
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Figure 49. The normalized arc length of the front-projected 
1/4-chord line of the flapping crane, seagull and goose wings 
as a function of time. 


Figure 47. The polynomial coefficients of the front-projected 
1/4-chord line of the flapping crane wing as a function of 
time. 
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Figure 50. The profiles of the front-projected 1/4-chord line 
of the flapping seagull wing at different instants. 



Figure 5 1 . The polynomial coefficients of the front-projected 
1 /4-chord line of the flapping seagull wing as a function of 
time. 



Figure 52. The orthographically projected semi-span b/2 
normalized by max(b/2) for the flapping seagull wing as a 
function of time. 



Figure 53. The profiles of the front-projected 1/4-chord line 
of the flapping goose wing at different instants. 
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Figure 54. The polynomial coefficients of the front-projected 
1 /4-chord line oi the flapping goose wing as a function of 
time. 
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Figure 56. An X-ray image of a seagull wing. 
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Figure 57. Two-jointed arm system. 
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Figure 55. The orthographically projected semi-span b/2 
normalized by ma x(b/2) for the flapping goose wing as a 
function of time. 
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rure 58. Projected views of 
) view, (b) side view, (c) the 
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Figure 62. Reconstructed flapping seagull wing at cot=0 . 
n / 4 . n / 2 , 3n / 4 and n . 



